Temperature Dependence of the Magnetic Susceptibility for Triangular-Lattice 
Antiferromagnets with spatially anisotropic exchange constants 
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We present the temperature dependence of the uniform susceptibility of spin-half quantum anti- 
ferromagnets on spatially anisotropic triangular-lattices, using high temperature series expansions. 
We consider a model with two exchange constants, Ji and J2 on a lattice that interpolates between 
the limits of a square-lattice (Ji = 0), a triangular-lattice (J2 = Ji), and decoupled linear chains 
(J2 = 0). In all cases, the susceptibility which has a Curie- Weiss behavior at high temperatures, rolls 
over and begins to decrease below a peak temperature, T p . Scaling the exchange constants to get the 
same peak temperature, shows that the susceptibilities for the square-lattice and linear chain limits 
have similar magnitudes near the peak. Maximum deviation arises near the triangular- lattice limit, 
where frustration leads to much smaller susceptibility and with a flatter temperature dependence. 
We compare our results to the inorganic materials CS2CUCI4 and Cs2CuBr4 and to a number of or- 
ganic molecular crystals. We find that the former (CS2CUCI4 and Cs2CuBr4) are weakly frustrated 
and their exchange parameters determined through the temperature dependence of the suscepti- 
bility are in agreement with neutron-scattering measurements. In contrast, the organic materials 
considered are strongly frustrated with exchange parameters near the isotropic triangular-lattice 
limit. 

PACS numbers: 75.10.Jm 



I. INTRODUCTION 



Understanding the interplay of quantum and ther- 
mal fluctuations and geometrical frustration in low- 
dimensional quantum antiferromagnets is a consider- 
able theoretical challenge 1^*2^1^ Research in frustrated 
quantum antiferromagnets was greatly stimulated by An- 
derson's "resonating valence bond" (RVB) paper— in 
which he suggested that the parent insulators of the 
cuprate superconductors might have spin liquid ground 
states and excitations with fractional quantum numbers, 
motivated by his earlier suggestion of such a ground 
state for the Heisenberg antiferromagnet on the trian- 
gular lattice— The Ising model on a triangular lattice 
illustrates the rich physics that can arise due to frustra- 
tion: it is known to have a macroscopic number of de- 
generate ground states— The antiferromagnetic Heisen- 
berg model with spatially anisotropic exchange interac- 
tions on the triangular lattice is of interest both theoreti- 
cally and experimentally. It describes the spin excitations 
in Cs2CuCL|i£ and Cs2CuBr^ii and the Mott insulating 
phase of several classes of superconducting organic molec- 
ular crystals^. Other materials for which this model is 
relevant include NaTiC^r 3, CuCl2 graphite intercalation 
compounds^ and the anhydrous alum, KTi(SC"4)2- 15 
Theoretically, this Heisenberg model is a candidate for 
a system with spin liquid ground states and possibly ex- 



citations with fractional quantum numbers— For the 
triangular-lattice model with spatially isotropic interac- 
tions, the preponderence of numerical evidencei&i^a^fliSI 
suggests that the ground state has long-range mag- 
netic order. However, making the interactions spatially 
anisotropic can lead to a very rich ground-state phase 
diagram 2 — 

The spatially anisotropic model, defined by a nearest- 
neighbor exchange constant J\ along one axis and J2 
along all other axes [see Fig. ^ , interpolates between the 
limits of square-lattice (Ji = 0), triangular-lattice (J2 = 
Ji) and decoupled linear chain (J2 = 0) limits. 23,24 It has 
been studied by spin wave theory^S 5 . series expansions^ 
large-N techniques^ slave fermionsj2£ Schwinger bosons 
with gaussian fluctuations^ and variational quantum 
Monte Carlo— - Quantum fluctuations are largest for 
Ji ~ O.8J2 and for Ji > 4J 2j 22 i 25 and so for these pa- 
rameter regions one is most likely to observe quantum 
disordered phases. 

From an experimental point of view, it is highly de- 
sirable to have a definitive way to determine the values 
of the exchange parameters for individual material sys- 
tems. Recently, it has been shown how for materials with 
relatively small values for the Heisenberg exchange con- 
stants J this can be done in high magnetic fields, using 
inelastic neutron scattering to measure the spin wave dis- 
persion in the field induced ferromagnetic phase— - The 
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temperature dependence of the magnetic susceptibility 
is one of the most common experimental measurements 
and it would be very useful if that can be used to deter- 
mine the extent of frustration and the various exchange 
constants directly. It is particularly important to have 
a scheme for materials, where the very high tempera- 
ture behavior of the system (T ^> J) is not accesible 
to experiments. Previously, Castilla, Chakravarty, and 
Emery pointed out how the temperature dependence of 
the magnetic susceptibility of the antiferromagnetic spin 
chain compound CuGeC"3 implied significant magnetic 
frustration^ In that case, it constrains the ratio of the 
nearest and next-nearest neighbour exchanges along the 
chains Similarly, it is reasonable to expect that the 
temperature dependence of the magnetic susceptibility 
should depend on frustration in two-dimensional models 
also and hence constrain the ratio Ji/J 2 . 

The Mott insulating phase of the organic molecular 
crystals is of particular interest because under pressure 
the materials considered become superconducting. A 
possible RVB theory of superconductivity in such mate- 
rials, emphasising the role of frustration, have recently 
been proposed^ These materials have exchange con- 
stants in the range of several hundred Kelvin, and their 
behaviour has led to several puzzles. Tamura and Kate 34 
measured the temperature dependence of the magnetic 
susceptibility for five organic molecular crystals in the 
family, /3'-[Pd(dmit)2]X (where dmit is the electron ac- 
ceptor molecule thiol-2-thione-4, 5-dithiolate, C3S5) and 
the cation X = Me 4 As, Me 4 P, Me 4 Sb, Et 2 Me 2 P, and 
Et 2 Me 2 Sb, where Me = CH 3 and Et = C 2 H 5 , denote 
methyl and ethyl groups, respectively). They compared 
their results with the predictions for the square and tri- 
angular lattices and found that for all the materials the 
results could be fitted by the high temperature series ex- 
pansion for the triangular lattice. However, some and 
not all of them undergo a transition to a magnetically 
ordered state at low temperatures. 

Recently, Shimizu et al.^ showed using *H nuclear 
magnetic resonance that k-(BEDT-TTF) 2 Cu 2 (CN) 3 did 
not undergo magnetic ordering and that the temperature 
dependence of the uniform magnetic susceptibility could 
still be fit by that for the triangular lattice. However, 
it should be stressed that for these molecular crystals 
the underlying triangular lattice of molecular dimers (to 
which each spin is associated) is not isotropic^ and so it 
is important to know the extent of the spatial anisotropy 
because this has a significant effect on the possible ground 
state. The isotropic triangular lattice is believed to be 
ordered, but for Ji/J 2 = 0.7-0.9 the anisotropic lattice 
could be quantum disordered^ Hence, determination of 
the actual ratio is important for understanding these ma- 
terials. 

Here, we use high temperature series expansions to 
calculate the temperature dependent uniform suscepti- 
bility of the spatially anisotropic triangular-lattice mod- 
els. Such calculations have been done previously for the 
pure square and triangular-lattice cases 3 ^! but not for 



the spatially anisotropic triangular-lattice model. This 
method is particularly useful here, as it allows one to 
cover the full range of Ji/J 2 ratios at once. Our main 
finding is that the susceptibility, for these antiferromag- 
nets, shows a broad maximum at a temperature (which 
we call the peak temperature T p ) of order the Curie- 
Weiss temperature. If the exchange constants are scaled 
to give the same peak location, the magnitude of the 
peak susceptibility varies with frustration. The unfrus- 
trated models, represented by the square-lattice and the 
linear-chain limits have similar peak susceptibilities. The 
triangular-lattice deviates the most from them, having a 
much smaller peak value, and a much flatter temper- 
ature dependence. The parameter regimes, where the 
ground states could be spin-disordered do not stand out 
in these calculations^ and are similar to the triangular- 
lattice limit. The reason for this is probably that at the 
temperature scales considered the susceptibility is largely 
determined by short-range frustration, rather than long 
length scale physics such as the existence of spin liquid 
states at zero temperature. 

Comparison with the measured susceptibility of 
Cs 2 CuCl4 and Cs 2 CuBr4 leads to exchange parameters 
in agreement with previous neutron measurements. For 
the organic materials, it shows that they are all close 
to the isotropic triangular-lattice limit. But, some of 
them could be weakly anisotropic, leading to a quantum- 
disordered ground state. Since, the organic materials are 
close to a Mott metal-insulator transition, we consider 
the possible role of multiple-spin exchange. Such inter- 
actions can be necessary for a quantitative description of 
such materials^ 



II. THE FRUSTRATED MODEL 

The spatially anisotropic triangular-lattice is shown in 
Fig. ^ The antiferromagnetic Heisenberg model is de- 
scribed by the Hamiltonian 

h = j x J2 s i • s ; + J 2 J2 s * ■ ^ w 

a b 

where the first sum runs over all nearest neighbor pairs 
along the x-axis and the second sum runs over all other 
nearest neighbor pairs. The vectors S represent spin- 
half operators. It is evident that, for J\ = 0, the model 
is equivalent to the square-lattice Heisenberg model, for 
J 2 = Ji it is equivalent to the isotropic triangular lattice 
model, and in the limit J 2 — > 0, it is equivalent to a model 
of decoupled linear chains. 

We now discuss how we might quantify how the 
amount of frustration in the model varies with </ 2 /Ji. 
Possible measures of frustration which have been dis- 
cussed before include: 

1. The number of degenerate ground states. 

2. How the competing interactions prevent the pair- 
wise collinear alignment of spins that would give neigh- 
boring spins the lowest interaction energy. 
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In order to quantify 2., Lacorre^ 9 . considered classi- 
cal spins and introduced a "constraint" function F c — 
—Eo/Eb which is the ratio of the ground state energy 
Eq of the system to the "base energy", Eb which is the 
sum of all bond energies if they are independently fully 
satisfied, i.e. 

Eb = — /J \Jij\(Si ■ Sj) max (2) 

Lacorre suggested that F c has values ranging from -1 
(no frustration) to +1 (complete frustration). However, 
for spin models that have a traceless Hamiltonian the 
ground state energy cannot become positive. So, F c must 
lie between -1 (unfrustrated) and (fully frustrated)- the 
largest possible value of F c is zero. Considering a single 
isosceles triangular plaquette taken from the lattice in 
Fig. ^ Lacorre found that for classical (large-.?) spins as 
a function of J2/J1, F c had its maximum value (—1/2) 
for the isotropic triangle (Ji = J2). The same result 
holds for the infinite lattice. 

Kahn 4 - recently stressed that for Heisenberg spins the 
degeneracy of the ground state depends on the value of 
the spin quantum number, 5, as well as the geometry of 
the plaquette. For example, on an isotropic triangle, the 
ground state is four- fold degenerate for 5 = 1/2 but non- 
degenerate for 5 — l4i On a single isosceles triangle, for 
5 = 1/2, the ground state has total spin 5r = 1/2 and is 
2-fold degenerate for J\ ^ J 2 and 4-fold degenerate at the 
isotropic point J\ — J 2 ■ We find that both F c is maximal 
(—1/3) and the ground state has the highest degeneracy 
for Ji = J 2 . On the other hand, for spin 5=1 the 
ground state is a non-degenerate singlet (St = 0) for 
a wide region near the isotropic limit (0.5 < Ji/J% < 
2), is three-fold degenerate (St = 1) outside this range 
(J1/J2 < 0.5 or J1/J2 > 2) and has accidental 4- fold 
degeneracy at the special points J1/J2 = 0.5,2. The 
function F c has no singular maximum, but a plateau at 
—0.5 for the whole range 0.5 < J1/J2 < 2, so the spin-1 
case is much less frustrated than the extreme spin-1/2 
case. 

The above properties of the degeneracy and constraint 
function are not unique to quantum spins but also hold 
for the Ising model on the same lattice. For a single 
isosceles triangle and for 5 = 1/2 the ground state energy 
changes at J\ = J 2 from — Ji/4 for J\ > J 2 , to (— 2J 2 + 
Ji)/4 for J 1 < J 2 . The "base energy" is E b = -(J\ + 
2J 2 )/A and hence F c has its maximum value (-1/3) when 
Ji = J 2 . The degeneracy of the ground state is 2 (only 
up-down symmetry) for J\ < J 2l 6 (only all up and all 
down are not ground states) for J± — J 2 , and 4 (either 
one of the J 2 bonds can be dissatisfied) for J\ > J 2 . So 
indeed by both measures for J\ — J 2 the model on a 
triangle is most frustrated. Extending this analysis for a 
single triangle to a large lattice of N sites the difference 
is even more dramatic as the degeneracy is£ exp(c-ZV) for 
Ji = J 2 , and is easily seen to be only 2 for J\ < J 2 
and exp(c'iV 1 / 2 ) for J\ > J 2 , where c, c' are numbers 



of order one. So the model has the largest ground-state 
degeneracy at the isotropic point. 

Although, this paper is concerned with the quantum 
spin model, the reason we mention the above properties 
of classical models is because an important question is 
whether our results concerning the connection between 
the amount of frustration and the temperature depen- 
dence of the susceptibility are also exhibited by the cor- 
responding classical Heisenberg and Ising models. This 
may be the case if the temperature dependence of the 
susceptibility down to the peak is largely determined by 
the frustration and correlations associated with a single 
placquette. 

With regard to measures of frustration we also note 
that from an experimental point of view two measures 
that have been proposed previously^ (i) The ratio of the 
Curie- Weiss temperature to the magnetic ordering tem- 
perature. This increases with increasing frustration, (ii) 
The amount of entropy at temperature scales much less 
than the exchange energy. 




FIG. 1: The spatially anisotropic exchange constants for the 
Heisenberg model on the triangular lattice. The model can 
also be viewed as a square lattice with an extra exchange 
along one diagonal. 



III. HIGH TEMPERATURE SERIES 
EXPANSIONS 

The high temperature series expansion method has 
been extensively applied to and tested for quantum lat- 
tice models^ We have obtained high temperature ex- 
pansions for arbitrary ratio of J\j ' J 2 to order (3 W . We 
express the uniform susceptibility, per mole, as 



where Na is Avogadro's number, g the g-factor, /is a 
Bohr-magneton, k the Boltzmann constant and T the 
absolute temperature. The dimensionless quantity \ can 
be expressed in a high-temperature expansion in J 2 /T 
and y — J\/ J 2 , as 

n 

X = E( J 2/ T )" E c m , n y m /(4 n+1 n!) (4) 

n— m — 



4 



The integer coefficients c m . n 
are presented in Table [I] 



complete to order n = 10 



IV. CURIE- WEISS BEHAVIOR AND BEYOND: 
SERIES EXTRAPOLATIONS 

As is well known, the high temperature behavior of the 
susceptibility, per mole, is given by a Curie- Weiss law 



C 



X 



T + T cw 

For our model, the Curie constant 

C = N A g 2 n%/4k = Ag 2 , 



(5) 



(6) 



with A — 0.0938 in cgs units. The Curie- Weiss tempera- 
ture is 



(7) 



From an experimental point of view, an important ques- 
tion is: How low in temperature is the Curie- Weiss law 
valid? To investigate this, we plot in Fig. the nor- 
malized inverse susceptibility as a function of T/T cw for 
several parameters, together with the Curie- Weiss law. 
It is clear that below T < 10T CW , the Curie- Weiss fit is 
no longer accurate. Deviations from the Curie- Weiss be- 
havior are the smallest near the triangular-lattice limit, 
and largest for linear chains. If one were to fit the inverse 
susceptibility below some temperature to a Curie- Weiss 
behavior, one would get a systematically larger Curie- 
Weiss temperature. To quantify this, we define an ef- 
fective temperature dependent Curie- Weiss constant T ( 
as 



eff 
cw 



t: 



eff 



X 



dx/dT 



(8) 



If one was to fit x _1 to a linear curve in the vicinity 
of some temperature (T) and use the intercept to esti- 
mate the Curie- Weiss constant, one would get T°%. Fig. 
Et" shows how T c e 5 varies with temperature for several 
parameter ratios. It shows that attempts to fit to a 
Curie- Weiss behavior below four times the Curie- Weiss 
temperature can result in an overestimate in the Curie- 
Weiss constant by less than 20 percent for the isotropic 
triangular lattice, whereas for the square lattice an error 
of 40 percent is possible. Similar observations were made 
previously for the classical Heisenberg antiferromagnet 
on a kagome latticed 

To obtain the susceptibility for T < T cw , we need to 
develop a series extrapolation scheme. We have used d- 
log Pade and the integral differential approximants to 
extrapolate the seriesMi2£i2&. For the linear chain model 
we use very long series given by Takahashi 4 ^ and for 
the square and triangular-lattice cases we have also used 
longer seriesi^S In the former case, the calculated sus- 
ceptibility agrees well with the exact results obtained 



u 





1 1 1 1 1 1 


1 1 1 1 | 1 1 




•— ■ Linear-Chain 






0—0 Square-Lattice 






»— « Triangular-Lattice 




'- 


<— > Jj/J 2 =3 






— Curie-Weiss Law 




1 , 


, I , , , 1 


1 ■ , , 1 , , ■ ' 



T/T 



i 

H° 1-4 



Linear Chain 
Square-Lattice 
Triangular Lattice 
J/J 2 =3 




T/T 



FIG. 2: (Color online) (a) Inverse magnetic susceptibility 
as a function of temperature in relative units of the Curie- 
Weiss constant T C w • The susceptibility departs from the high- 
temperature Curie- Weiss limit (Eq.(4)) already at tempera- 
tures a few times T cw due to short-range correlations. The 
smallest departure occurs for the triangular lattice, (b) Ef- 
fective Curie- Weiss constant vs temperature found by a 
local fit of the susceptibility to the Curie- Weiss form, Eq.(7). 
The plot shows that fitting data below 4T CW can result in large 
overestimates of the Curie- Weiss constant. 



from Thermodynamic Bethe Ansatz calculations^. For 
the square and triangular-lattice cases it also agrees well 
with previous numerical calculations- 6,4 — . In all cases, 
several integral/dlog-Pade approximants are calculated, 
and in the plots below two outer approximants are shown, 
i.e. a large number of approximants lie between those 
shown. Based on our general experience with series 
extrapolations^, we feel confident that as long as the 
upper and lower curves are not too far from each other, 
they bracket the true value of the thermodynamic suscep- 
tibility. In general, we find that the extrapolations work 
well down to the peak temperature and begin to deviate 
from each other below the peak. It is not possible to ad- 
dress the zero and very low-temperature behavior of the 
susceptibility from these calculations. 

In Fig. |3| we show the uniform susceptibility, for dif- 
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FIG. 3: (Color online) Susceptibility versus temperature for 
different values of J1/J2. The peak susceptibility \p an d the 
Curie- Weiss constant C = Ag 2 are used to define a dimension- 
less relative temperature scale. As discussed in the text, the 
two curves shown for each J1/J2 value are due to different 
extrapolation schemes. For the most frustrated triangular- 
lattice the peak in the susceptibility occurs at the lowest rel- 
ative temperature. 



ferent y = J1/J2, as a function of temperature. For 
all J1/J2 ratios, there are two plots showing the upper 
and lower limits of extrapolated values as discussed in 
the previous paragraph. The susceptibility is scaled to 
have a peak value of unity, and the temperature axis is 
scaled by the peak susceptibility to a dimensionless rela- 
tive temperature. One finds that the susceptibility peaks 
at a comparable relative temperature for the unfrustrated 
square-lattice and linear chains. The primary difference 
between these two models lies in the behavior of the sus- 
ceptibility below the peak. It decreases much more slowly 
for the linear chains than it does for the square-lattice. 
We believe, this is related to the fact that longer-range 
antiferromagnetic correlations grow much faster for the 
square-lattice than they do for linear-chains. Thus the 
shift of the spectral weight away from zero wavevector oc- 
curs more gradually for linear chains. For the triangular- 
lattice, the peak is shifted to much lower relative tem- 
peratures. Note that the triangular-lattice has a peak 
at a temperature even lower than for J1/J2 = 0.8, where 
T = calculations show an absence of long range order—. 

From Fig. [21 it is clear that frustration leads to a re- 
duction in the magnitude of the product x p T p as well 
as a reduction in the peak temperature T p with respect 
the Curie- Weiss temperature T cw . These parameters are 
plotted in appropriate dimensionless units in Fig. 01 as 
a function of the frustration ratio J\/(J\ + J2), and 
both have a minimum around the triangular lattice limit 
J 1 = J 2- To connect with experiments we also show the 
ratio of the peak temperature T p and the Zeeman energy 
required to fully polarize the spins <?/XB-B sa t, related to 




0.4 0.6 
J 1 /(J 1 +J a ) 

FIG. 4: (Color online) Variation of the key parameters of the 
susceptibility curve as a function of the frustration ratio 
Ji/(Ji + J2): location of peak temperature T p relative to the 
overall energy scale of the couplings, given by the Curie- Weiss 
constant T cw or the saturation field B sat required to overcome 
all antiferromagnetic interactions (see text for more details), 
dimensionless product of peak susceptibility and peak temper- 
ature TpXp/Ag 2 (with A — .0938 in cgs units), and flatness 
of the susceptibility curve x(4Tp) / 'x(T p ) . The isotropic trian- 
gular lattice (Ji = J2) is the most frustrated with the lowest 
relative peak temperature T p /T c „, lowest peak susceptibility 
and flattest curve at temperatures above the peak. The circles 
are values extracted from the experimental data for CS2CUCI4 
from Ref. [m] and the squares are for Cs 2 CuBr 4 Ref. El This 
suggests that the ratio J1/J2 is close to 3.0 and 2.0, respec- 
tively, for these two materials. The former is consistent with 
independent estimates from neutron scattering^ 



the couplings strength 




•'2 

2Ji 



for J 2 < 2Ji 
for J 2 > 2Jx . 



(9) 



Fig. 21 also shows the ratio x(4r p )/x P , which is a mea- 
sure of the flatness of the curves on the high-temperature 
side of the peak. A larger value of this ratio implies 
a slower decay of the susceptibility with temperature. 
These quantities clearly show that the triangular-lattice 
is the most frustrated, with the lowest peak temperature 
relative to the scale of the exchange interactions, T p /T cw 
or kTp/gfiBB sa t, the smallest dimensionless ratio T P \ P 
and the flattest peak denoted by the largest x(4T p )/x P - 
The plots look very symmetrical around the triangular- 
lattice limit, and there is nothing anomalous about the 
case of J1/J2 = 0.8, where zero-temperature studies give 
a disordered and gapped dimerized ground stated. We 
note that all of the extracted parameters in Fig. 0] are 
from the susceptibility curve at temperatures above the 
peak and in order to see evidence for the presence of a 
gap for J1/J2 ~ 0.8 as opposed to no gap in the isotropic 
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triangular-lattice case one would be required to analyze 
the susceptibility curve at temperatures much below the 
estimated gap A ~ 0.25J2 ~ 0.5 T p in the dimerised 
stated and such low temperatures are not accessible by 
the present series calculations. 



V. COMPARISON WITH EXPERIMENTAL 
SYSTEMS 

In this section, we compare our theoretical results with 
experimental data on CS2CUCI4, Cs2CuBr4 and various 
organic materials. In Fig. [3 we show the susceptibility as 
a function of temperature for different J1/J2 ratio, where 
the temperature is scaled by the peak temperature (T p ) 
and the susceptibility itself is scaled by the peak tempera- 
ture to give a dimensionless reduced susceptibility. This 
plot is very instructive as it allows one to clearly read 
out the J1/J2 ratios. Also shown are the susceptibilities 
for the materials CS2CUCI4 and Cs2CuBr4, with their 
g-values taken from ESR experiments^^^. In this plot 
with no free parameters, it is apparent that the J1/J2 ra- 
tio is near 3.0 for CS2CUCI4 and near 2.0 for Cs2CuBr4. 
Some of these results can also be seen from Fig.0J where 
key dimensionless ratios of the temperature dependent 
susceptibility are shown. 

A more detailed comparison of the susceptibility for 
the materials, CS2CUCI4 and Cs2CuBr4, allowing g to 
vary freely is shown in Fig. Once g is allowed to vary, 
the material CS2CUCI4 can be fit above the peak not 
too badly even with the pure square-lattice model (not 
shown). However, a much improved fit happens with 
J1/J2 = 3 and J 2 = 1.49 K in excellent agreement with 
the exchange values extracted directly from neutron scat- 
tering measurements^ Also shown are fits to linear chain 
and triangular-lattice limits, which bracket J1/J2 = 3. 
One can see significant deviation in both limits. The 
large deviation from the isotropic triangular-lattice case 
shows that frustration is relatively weak in this material. 

For Cs2CuBr4, the best fit for J1/J2 = 2 arises with 
J 2 = 6.99 K. However, when g is allowed to vary, a 
range of J1/J2 values from 1.8 to 2.8 give comparable 
fits, several of which are shown in figure. In general, 
the high temperature data is better fit by a larger Ji / J 2 
value, whereas the data at and around the peak is bet- 
ter fit by a smaller J1/J2 value. No choice of parame- 
ters can fit the very low temperature data (below half 
the peak temperature). These values are also consis- 
tent with previous estimates. Using the value of the in- 
commensurate ordering wavevector Q = 0.575(l)b* ob- 
served by neutron scattering^ classical spin- wave theory 
gives J1/J2 = 2.14 whereas including quantum renor- 
malization corrections as predicted by large-A Sp(A) 
theory^ gives J1/J2 ~ 1-8, and series expansions^ gives 
J1/J2 ~ 1-4. This calls into question the rather large 
renormalization of the ordering wavevector found in the 
series expansion study. 

Now we turn to the organic materials. In Fig. {7\ we 
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FIG. 5: (Color online) Susceptibility vs temperature in units 
of the peak temperature T v . The isotropic triangular lattice 
(green line) has the lowest and flattest susceptibility. Solid 
squares show data points for the anisotropic triangular lattice 
material CS2CUCI4 (a- axis, Ref . IHll) . solid circles show data 
for Cs 2 CuBr 4 from Ref. El 
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FIG. 6: (Color online) Fits of the susceptibility in CS2CUCI4 
(a) and Cs2CuBr4 (b), see text. 
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show a corresponding comparison for the material k- 
(BEDT-TTF) 2 Cu 2 (CN) 3 . Only the theoretical data for 
the isotropic triangular-lattice are shown. One can see 
an important difficulty in using the T p -scaled plots near 
the triangular-lattice limit to determine J1/J2. For the 
organic material, k-(BEDT-TTF) 2 Cu 2 (CN) 3 , the mea- 
sured susceptibility is very flat and it is difficult to de- 
termine the peak temperature T p . From the data, the 
peak temperature appears to be between 65 K and 95 K. 
Using the values for T p of 65 K and 95 K, one can either 
get the data to fall above or below the triangular-lattice 
values. A suitably chosen peak temperature allows one 
to get very close agreement with the triangular-lattice 
limit. This peak temperature can also be used to deter- 
mine the exchange constant. However, for the triangular- 
lattice, there is theoretical uncertainty in the peak loca- 
tion. Hence, it is more accurate to directly fit the ex- 
perimental data to theory to obtain the exchange con- 
stants. For k-(BEDT-TTF) 2 Cu 2 (CN) 3 , fixing g = 2.006 
and Ji = J 2 , the best fit leads to J\ = 256 K, a value 
close to that obtained by Shimizu et alm^ 

The ability to fit flat susceptibilities to the isotropic 
triangular-lattice model is further illustrated in Fig. [S] 
where the susceptibility data are shown from five dif- 
ferent molecular crystals in the family /3'-[Pd(dmit) 2 ]X 
(where dmit is the electron acceptor molecule thiol-2- 
thione-4, 5-dithiolate, C3S5) and the cation X = Me4As, 
Me 4 P, Me 4 Sb, Et 2 Me 2 P, and Et 2 Me 2 Sb, where Me = 
CH3 and Et = C 2 Hs, denote methyl and ethyl groups, re- 
spectively). We have taken the g- value to be 2.04. By ad- 
justing the peak temperature, they can all be brought to 
rough agreement with the triangular-lattice model. As- 
suming isotropic interactions, and g = 2.04, we estimate 
the exchange constants to be 283 K, 289 K, 270 K, 279 K 
and 247 K respectively. It is clear that none of these 
organic materials are far from the isotropic triangular- 
lattice limit. But, we emphasize, that by this method 
it is difficult to discriminate between J\ / J 2 ratios in the 
range .85 < Ji/J 2 < 1.15. Note that the latter regions 
also include quantum disordered phases. 

To avoid the problem of determining the peak temper- 
ature, we go back to Fig.|3 and scale the data by the peak 
susceptibility. These can be inferred accurately from the 
data, even when the peak temperature cannot. In Fig. 

we show such a comparison of experimental data with 
theory. The data for k-(BEDT-TTF) 2 Cu 2 (CN) 3 lies ex- 
tremely close to the isotropic triangular-lattice case. The 
other materials deviate from the J\ = J 2 limit, but 
still lie in the range 0.85 < Ji/J 2 < 1.15. If we as- 
sume that the systems are described by the isotropic 
triangular-lattice, the exchange constant can be read 
of from the peak susceptibility by using the relation 
J = 0.0035g 2 /xp- This leads to exchange constants of 
250K for k-(BEDT-TTF) 2 Cu 2 (CN) 3 and 280 K, 289 K, 
260 K, 273 K and 236 K for the other materials. These 
values are close to those obtained from the best fits. 

It should be noted here that in the experimental data, a 
Curie term from magnetic impurities and a diamagnetic 



term has been subtracted and these can also influence 
the determination of exchange parameters. However, it 
is unlikely that any of these materials are very far from 
the isotropic triangular lattice limit. 

From the fits the Heisenberg couplings are compara- 
ble for all materials and around 250 K. We now consider 
how these compare with quantum chemistry calculations. 
The exchange constants can be related to parameters in 
an underlying Hubbard modeiiSiSi^i where J = 2t 2 /U 
and t is the intersite (i.e., inter-dimer) hopping and U is 
the cost of double occupancy for two electrons or holes 
on a dimer. If the Coulomb repulsion Uq on a single 
molecule within the dimer is much larger than the inter- 
molecular hopping to within a dimer then U ~ 2io- For 
/3'-[Pd(dmit) 2 ]X electronic structure calculations based 
on the local-density approximation (LDA)^^^ give 
t ~ 30 meV and t ~ 500 meV, and so J - 50 K. For k- 
(BEDT-TTF) 2 Cu 2 (CN) 3 Huckel electronic structure cal- 
culations give t ~ 50 meV and t ~ 200 meVr^. The 
resulting U ~ 400 meV is comparable to that deduced 
from measurements of the frequency dependent optical 
conductivity of similar k materials ^2*S This value of U 
is smaller than values deduced from quantum chemistry 
calculations on isolated dimers, which do not take into 
account screening^ Using the above values of t and U 
gives J ~ 100 K. 




^k-(BEDT-TTF) 2 Cu 2 (CN) 3 , T p =65K 

b-h Triangular-Lattice 



FIG. 7: (Color online) Comparison of the temperature 
dependence of the magnetic susceptibility of k-(BEDT- 
TTF)2Cu2(CN)3 with series expansions calculations for 
isotropic triangular-lattice. The experimental data is from 
Ref. |23 A value of g=2.006 was used based on electronic 
spin resonance measurements^ We see that this material is 
well described by a Heisenberg model on the isotropic trian- 
gular lattice, with peak temperature T v — 85 K. Note also 
that that the agreement is quite sensitive to changes in the 
value of T p , a quantity that is difficult to pin-point in a flat 
curve. 

Note that the quality of fit is best for the material n- 
(BEDT-TTF) 2 Cu 2 (CN) 3 , where it really fits well with 
the isotropic triangular-lattice model. However, it is 
also a system that does not order down to very low 
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FIG. 8: (Color online) Comparison of the temperature depen- 
dence of the magnetic susceptibility of five different organic 
molecular crystals from the family /3'-[Pd(dmit)2]X (different 
X are indicated with Me=CH3, Et=C2Hs) with series expan- 
sions for isotropic triangular-lattice. Experimental data is 
from Ref. I34I A value of g = 2.04 was used based on elec- 
tronic spin resonance measurements iS- All of these materials 
are well described by a Heisenberg model close to that for 
the isotropic triangular lattice, assuming that ring-exchange 
interactions do not need to be taken into account. 




TX D /Ag 



then these terms may be important. Recent neutron scat- 
tering studies showed the effect of such interactions on 
the dispersion of spin excitations in La2Cu04^ 

The metallic phase of the organics are in the regime 
U/t ~ 5 — 10 [Ref. 0| and so one might expect multiple- 
exchange terms to be relevant in the insulating phase. 
For the triangular lattice triple exchange is also pos- 
sible. However, for spin-1/2 this just corresponds to 
a renormalisation of the nearest-neighbour two-particle 
exchanged The frustrating effects of multiple-spin ex- 
change on the isotropic triangular lattice lead to rich 
physics and have an experimental realisation in monolay- 
ers of solid 3 He on graphite^ Let J denote the nearest- 
neighbour exchange and J4 the multiple spin-exchange, 
involving the four spins comprising a pair of triangular 
placquettes. This model has been studied extensively and 
exact diagonalisation calculations suggest that the 120 
degree Neel state, which is the ground state for J4 = 0, 
is destroyed when J4 > O.lJp^ It is appealing to think 
that this could be the explanation for why k-(BEDT- 
TTF)2Cu2(CN)3 does not magnetically order, whereas it 
should if it is really described by the isotropic triangular 
lattice nearest neighbour model. This material is close 
to a Mott-Hubbard metal-insulator transition since the 
insulating state is destroyed under pressure or uniaxial 
stress^SiSi However, it is not clear that J4 will be large 
enough in the actual material. The expressions derived 
from a t/U expansion giveSi J 4 /J = 10(t/t/) 2 . This 
means one must have U < St to obtain a spin liquid. 
However, exact diagonalisation of the Hubbard model on 
the isotropic triangular lattice at half filling, shows that 
the insulating state only occurs for U > 12t~ Hence, it 
is not clear that multiple-spin exchange could account 
for the fact that this material appears to be close to 
the isotropic triangle but does not magnetically order. 
However, to definitively resolve this issue would require 
a detailed study of the spatially anisotropic model with 
four-spin exchange. 



VI. CONCLUSIONS 



FIG. 9: (Color online) Parameter free comparison of the 
susceptibility data on organic materials with theoretical plots 
scaled by the peak susceptibility, Xp> which is easy to measure 
accurately for a flat curve. It is evident that the materials de- 
viate only slightly from the isotropic triangular-lattice model 
and have J1/J2 ratios in the range 0.85 to 1.15. 



temperatures^ This remains a puzzle. The quality of 
fits was not as good for the other organic compounds. 
It is quite possible that the organics have other interac- 
tions not captured by the Heisenberg model. In a Mott 
insulator when a perturbation expansion in t/U is used 
to derive an effective Hamiltonian for the spin degrees 
of freedom one finds that to four order in t/U there are 
cyclic exchange terms in the Hamiltonian iSi If U/t < 10 



In this paper, we have developed high temperature ex- 
pansions for the uniform susceptibility of the spatially 
anisotropic triangular lattice Heisenberg model. We find 
that the temperature dependence of the susceptibility at 
temperatures of order the exchange constants, are sen- 
sitive to frustration, that is, the ability of spins to align 
antiparallel to all their neighbors. The square-lattice and 
linear chain limits have similar reduced susceptibilities at 
and above the peak, while the triangular-lattice limit ap- 
pears most frustrated, with the smallest and flattest sus- 
ceptibilities. Comparison with various experimental sys- 
tems shows that a variety of organic materials are close 
to the isotropic triangular-lattice limit, whereas the inor- 
ganic materials CS2CUCI4 and Cs2CuBr4 are much less 
frustrated. 

It would be nice to have a simple formalism which 
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could provide an analytic relation between the peak sus- 
ceptibility and exchange parameters. Qualitatively, our 
arguments show that short-range frustration, or the in- 
ability to align parallel with respect to neighbors as quan- 
tified by the parameter F c in Section II, is maximum 
near the isotropic limit and this is what pushes the peak 
in the susceptibility down to lower temperatures. For 
a wide range of frustrated antiferromagnets it has been 
previously pointed out that the Curie- Weiss law holds to 
relatively low temperatures Several theoretical mod- 
els, mostly for classical spins, have been developed to ex- 
plain thisSiSi. Basically, frustration leads to individual 
plaquettes or spin clusters behaving essentially indepen- 
dently. However, our models are less frustrated than that 
and hence always develop substantial correlations. This 
means that any simplistic explanation is unlikely. 

In organic molecular crystals a weak temperature de- 
pendence of the magnetic susceptibility is often inter- 
preted as being evidence for metallic behavior, since for a 
Fermi liquid the susceptibility is weakly temperature de- 
pendent. However, this is inconsistent with the fact that 
in most of these materials above temperatures of about 50 
K there is no Drude peak in the optical conductivity and 
the resistivity has a non-monotonic temperature depen- 
dence and values of order the Mott limit?- 12-70 This work 
shows that due to the substantial magnetic frustration 
the susceptibility can actually be due to local magnetic 
moments, even though in the range up to 300 K one does 
not see a clear Curie temperature dependence. 

In a future study we will consider the temperature de- 



pendence of the specific heat capacity for this model. A 
previous study^i of the square lattice, single chain, and 
triangular lattice Heisenberg model found that the peak 
in the specific heat versus temperature curve occurred 
around J for all models but was much broader for the 
triangular lattice. A related issue was that as the temper- 
ature decreases the entropy decreases much more slowly 
for the triangular lattice than the others. 
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TABLE I: Series coefficients for the high-temperature expansions of the uniform susceptibility \ m Eq. (|1J. Nonzero coefficients 
up to order n = 10 are listed. 
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